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Abstract 

We present 0{a1) corrections to the decay t — > bW in the hmit of a very 
large top quark mass, mt ^ my/- We find that the 0{a^g) effects decrease the 
top quark decay width Tt by about 2%: Vt = To{l — 0.8as(mt) — l.lal). The 
complete corrections are smaller by about 24% than their estimate based on 
the BLM effects 0{l3Qa1). We explain how to compute a new type of diagrams 
which contribute to T{t —>■ bW) at the 0{a1) level. 



1 Introduction 



Because of its very large mass, the top quark is an interesting object for precise 
studies and searches for possible "new physics." In consequence of the large m^, top 
lifetime is very short. With good accuracy the top decay rate Tt is proportional to 
the third power of the top mass and equals approximately 1.5 GeV for rrit = 175 
GeV. This large width makes the top quark behave almost like a free quark, a unique 
situation in QCD ||T[|. The reason for this is that the time scale of the weak top 
quark decay 

[8V271 

is much shorter than the time scale of the non-perturbative QCD effects, 1/Aqcd- 

Within the Standard Model (SM) the reaction t bW is the dominant top quark 
decay channel. However, in most extensions of the standard model (SM) there 
are additional decay channels (e.g. t —>■ bH^ or decays into final states including 
super symmetric particles). Since the SM prediction for the top quark lifetime can 
be given with high precision, its measurements can in principle shed light on the 
new physics contributions to top decay. 

The fact that the life time of the top quark is so small has rather interesting conse- 
quences for the top production in e~^e~ annihilation. Because the nonperturbative 
effects are suppressed by the large Ff, threshold production cross section of ti in 
e~^e~ collisions can be predicted rather reliably. Clean signature and a rich physics 
potential make the investigation of this process in the threshold region one of 
the priorities at the Next Linear Collider (NLC). Many studies have recently been 
performed in order to describe this reaction with high precision (see e.g. @, |^). 

The total cross section for e'^e~ tt in the threshold region is determined by the 
equation: 

atis) r^lmG{E + iTt;0,0), (1) 

where s is the center of mass energy squared, E = ^/s — 2mt, and G{E] 0, 0) is the 
value at the origin of the Schrodinger equation Green function describing the non- 
relativistic ti system. Because of the (Tt{s) dependence on Tt, a precise prediction 
for the cross section requires a matching accuracy in Tf. This is true not only for 
the total cross section but also for other top threshold observables, for example top 
momentum distribution. 

The analysis of crt{s) has been recently extended to the next-to-next-to- leading order 
accuracy ^. That analysis was made possible in part by the recently calculated 
two- loop corrections to the QCD potential 0, whose influence on the threshold 
cross section was investigated in 0. On the other hand, the 0{al) corrections to 
the width were not known at that time and were not included into those studies. 
The purpose of the present paper is to provide Tt including C(af) effects. 

The direct measurement of the top quark decay width at NLC is difficult. At the 
moment a determination of the top quark width based on the measurement of the 
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forward-backward asymmetry of t quarks in the threshold region in e~^e~ colhsions 
is considered to be the most promising option. It is beheved that the accuracy of 
10 — 20% can be obtained in such a measurement p|. More optimistic estimates of 
5% accuracy have also been published 0. Still better accuracy might be obtained 
at a future muon collider. 

Because of the importance of the top quark width, much effort has been invested into 
studies of radiative corrections to its decays. Here we briefly summarize the results 
obtained in the SM. For a more extensive discussion and a summary of calculations 
in some extensions of the SM we refer the reader to ref. |]^, |TD|, where also corrections 
to various differential decay distributions are presented. 

QCD corrections to heavy quark decays were first studied in an effective Fermi-like 
theory, valid for quarks much lighter than the W [^. For the semileptonic decays 



such calculations were technically similar to the muon decay case |]T2|. After it had 
become clear that the top quark is heavy, QCD corrections were calculated in 
with full W propagator taken into account. These results are also applicable to the 
process t — ^ bW which is now considered the dominant decay channel of the top 
quark, and in this context were confirmed in []I4[. The one- loop QCD corrections 
decrease the rate of this decay by about —8.4% (for ~ 175 GeV). 

At the Olal) level only the so-called BLM [|15| corrections have been known |1l6| , |17|. 



Their numerical value will be illustrated below. Summation of these effects to all 
orders has been discussed in O, |l9ll. 



The electroweak corrections to t — bW were evaluated in ||20| and were found to 
increase the decay rate by about 1.7%. This effect is almost canceled by accounting 
for the finite width of the W, which decreases the rate by about —1.5% [l^. In the 
present paper these both effects will be neglected. 

Present uncertainty in the theoretical prediction of the top quark width is to a large 
extent connected with the unknown two-loop QCD corrections Q and the need for 
their complete evaluation has been repeatedly emphasized. This paper is devoted 
to this calculation. An exact computation of this effect would be very difficult, but 
since the size of those corrections is expected to be small, it is justified to make 
some approximations. In the present calculation we neglect the mass of the W; the 
error is expected to be of the order m^^/mf ^ 0.2.|^ Another, smaller source of error 
is inherent in the method we employ in this calculation. As described below, we 
use the difference of the t and b masses as an expansion parameter and calculate 
many (about 20) terms of the series in this parameter. An error of up to 6% is 
caused by the truncation of this series. In any case, the accuracy of our final result 
is completely sufficient for all applications which can presently be contemplated. 

In this paper we consider W as stable. Part of the effects of its instability can 
be accounted for using the known results for the 0{al) corrections to its hadronic 
width. There are also "non-factorizable" corrections due to an exchange of gluons 

^At the 0{as) level, the difference between the QCD correction calculated for mw — rnb = 
and the correction calculated for physical values of niw and mt is approximately 10%. 
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between the top, or the b quark resulting from its decay, and the quarks produced 
in the W decay. We expect these effects to be additionally suppressed by Tw/rriw 

There is also another application of the calculations reported in this paper, related 
to the b quark physics. The inclusive decay width b — >■ ului is of interest for the 
extraction of the CKM matrix parameter \Vub\- The 0(a;^) corrections to t — > bW 
presented below, after simple modifications can be used to obtain the corrections to 
the differential inclusive decay width r(6 —>■ ulvi) at the point where the invariant 
mass of leptons vanishes. Thus, the results reported here can be used in the future 
to estimate the second order QCD corrections to 6 ^ ulvi. 

Most of the results presented here have been derived using methods described in ||21| 
in the context of 6 — >■ c decays. In addition, since the t quark is much heavier than 
6, we have to include three diagrams which result from "non-planar" interference of 
t Wbbb amplitudes. Their calculation is described in some detail in the following 
section |^. In section ^ we present the numerical results for various values of the 
quark masses in the final state. Our conclusions are given in section ^. 

2 Non-planar diagrams t bbbW 

Since the top quark is significantly heavier than myy + Sm?,, the decay t bbbW 
contributes to the width of the top at 0{al). The treatment of this decay channel 
is complicated by the presence of non-planar interference diagrams shown in fig. ^ 




Figure 1: Non-planar diagrams contributing to the decay channel t — > Wbbb. 

In order to compute them, we expand around the point rrib = rrit/S. The expansion 
variable is 6 such that 

m, = ^(l-5). (2) 

We introduce the following notations: pi^2 the momenta of the virtual gluons which 
later split into bb pairs, namely: pi = + P5 and P2 = Ps + Pa- P3,P5 are the b 
quark momenta and p4 is the antiquark momentum in the final state. W denotes 
the VT-boson momentum, t is the top quark momentum, p^ and Pq are the momenta 
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of the virtual top propagators: Pe — t — pi and Pr = t — p2- b is the momentum of 
the virtual 6-quark, b = t — W . 

All propagators can be expanded around the static limit, because the configuration 
around which we expand is the top quark decaying into three 6-quarks at rest. For 
instance, the leading term in the expansion of the gluon propagator l/pl is l/(4m^). 
The virtual top propagators give: 

1 - ' (3) 

{t - piY - ml ^mb{mt - nib) ' 

etc. Therefore, no dependence on phase space variables remains in the expanded de- 
nominators and the integration over the phase space can be performed by successive 
factorization into two-particle phase spaces. 

First we introduce the vector H which combines p4 and p^: H — p^ + p^. We 

decompose ^4^5 into components parallel and perpendicular (in four-dimensional 
sense) to the H direction. After averaging over directions of the perpendicular 
components this phase space integration gives the volume factor 



R,(H;p,,p,) = -^l-^. (4) 

Next, we combine H and W into Q — H + W. Since W is masseless, this phase 
space gives no square root function 

R2{Q;H,W)^'^^^^. (5) 

We are then left with t — Q+ps phase space and two integrations over H"^ and Q^. 
The integration limits are 4m^ < Q'^ < {M — m)'^ and 4m^ < H"^ < Q"^. Also, the 
Q,P3 phase-space is given by Kallen function, the momentum of the Q in the rest 
frame of t: 

R2{t;Q,P3) = ^VfK - ^f)' - + ^^')' - Q']- (6) 



Finally, for and we introduce variables xi^2 G (0, 1), defined by = 4m^(l-|- 
10x1X2) and Q'^ = 4m^(l -|- uxi). From the limits on H'^ and Q'^ we find 



UJ — 



{m t - rribf 
Ami 



- 1 ~ 5 (7) 



and the volume of the phase space R2{t; Q,P3) is given by 



R2{t;Q,P3) = ^V^VT^y 1 - J - ^(4 - ~S), (8) 



where the last square root can be expanded in 5. Also the y 1 — Am'^/H'^ expressed 
in terms of 0:1^2 becomes proportional to y/XiX2- As a result we obtain an expression 
containing only half-integer or integer powers of Xi 2 which can be easily integrated. 
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3 Two-loop corrections to top decay rate 



In the limit of vanishing W mass, the width of the decay t 
formula (we neglect terms 0{a^)) 



bW is given by the 



T{t bW) = To 



ao H Cpai + 



TT 



(9) 



where Tr 



Gprnf 



Vfftp. In this paper we use as{mt) determined in the MS scheme 



and mt b are the pole masses, 



ml/mf)^ and ai is also known in a closed 



analytical form |]T3l In the limit of the vanishing b quark mass we have 



ao = 1, 



_ 5_7r2 



(10) 



a2 is the main new result of the present paper. It can be divided up into 4 gauge 
invariant parts 



a2 



Cf y) + ^AttA + TjiNiaL + Tjian, 



where Cp = 4/3, Ca = 3, Tr = 1/2, and Nl is the number of quark flavors whose 
masses can be neglected. For the top quark decay we take A''^ = 5. ap^A,L,H are 
functions of the mass ratio of the b and t quarks. We have calculated them using an 
expansion around the equal mass case. The expansion parameter, 6 = 1 — mh/rrit is 
close to one, so that many terms of the series must be calculated in order to obtain 
good accuracy. Such expansion has already been considered in our previous work 
23| , In those papers we have defined the quantities Af,a,l,h- Two coefficient 
functions in eq. (|l3), cla,l, coincide with (5^A^ j;^.^ For the purpose of the present 
application we have more than doubled the number of terms in the expansion, 
is analogous to ^^A^^, except that it now contains contributions of the top quark 
loop only. The b quark loop and part of real bb pair radiation is accounted for by 
adjusting Nl. There is, however, part of real bb pair production which cannot be 
described in this way. It corresponds to the non-planar interference diagrams shown 
in fig. and is included in ap, which without these diagrams would have coincided 
with S^Ap. These diagrams, contributing to the decay t Wbbb, have no analogs 
in the b — c transitions considered in [^, 

For nib 0, ol was obtained numerically in |T^, and is now known analytically in 
this limit 



17 : 



a limb = 0) 



4 

-9 + 



23 
T08' 



-IT 



+ C3 



Cs ^ 1.2020569. 



(12) 



In order to obtain the correct result one has 

2 /^si Y^oo fl„+4^ 



^ There is a mistake in the formula for A a in |l[ 
to add to 6^Aa given there the following expression: — ^(Cs — 7r^/6) X^^i f}"'~^'^n/[{n + l){n + 2)] 
with P = (5(2 — (5), expanded in 6 to desired order. For b — > clvi the resulting change is insignificant; 
the magnitude of Aa is decreased by less than 3%. 
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The dependence of the coefficient functions a, on mi,/ nit is shown in fig. |^. For 
rrib/mt ^ 4.8/175 we find the following numerical values: 

For rrib/mt = 4.8/175: 

ap = 3.4(2), ttA = -6.28(2), 

ttL = 2.83(4), an = -0.06349(1). (13) 

We note that ap has the largest error bar. This is because there are two separate 
series contributing to it: in addition to the series in 5 = 1 — m^/mt there is also 
the non-planar contribution expanded in S = 1 — 3mf,/mj. Both these series are 
separately divergent in the limit mf, — > like In^(mb). For this reason their reliable 
estimate near this limit is difficult and we have assigned a conservative error bar 
to the result. On the other hand, because of the color factors in eg . ([TT|) , the final 
result for 02 is not very sensitive to the relatively large error in a^. 

The expressions we have obtained for ap^A^i contain the ffist and second powers of 
ln(25). It turns out that the convergence of these series is improved if we rewrite 
these terms as \n{26) = ln/5 — ln(l — 5/2), with (3 = 6{2 — 6) and expand ln(l — 5/2) in 
6. Since [3 approaches 1 much faster than S, the remaining logs are smaller. Because 
of the many terms which had to be evaluated for the coefficient functions Oj, the 
results are rather lengthy and not suitable for publication in a journal. However, 
they can be obtained from the authors upon request. 

In the future, it is likely that the corrections to t bW will be calculated analyt- 
ically in the point = 0. We therefore list here our estimates for the numerical 
values of in this point: 

For nil, l^t = : 

ap = 3.5(2), ttA = -6.35(7), 

aL = 2.85(7), an = -0.06360(1). (14) 
The value of ai found from our expansion is in agreement with the exact value 



obtained from eq. ([ID, aLirub = 0) = 2.859473 , 



Let us now discuss the numerical value of the correction in the case of the massless 
b quark. For the purpose of discussion it is convenient to separate the BLM |TB[ and 



the non-BLM corrections to the decay rate. The BLM corrections follow immediately 



from the results of |]T6|, and in present notations are given by: 



^ = -T^, (f - iV,) a,. (15) 

Using Nl = 5, we obtain: 

afM^^^^O) = -16.442.... (16) 
This is to be compared with the complete result for 02, 

a^{mf, = 0) = -12.5(4). (17) 
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Figure 2: Plots of the coefficient functions defined in eq. (|n|): (a) ap, (b) a^, (c) 
ai, (d) an- The variable on the horizontal axis is nib/ nit. 



The complete 0{a1) corrections are smaller by 24% than the BLM estimate. This 
difference is somewhat larger than in 6 — > c decays p3|, ^ . 



In any case, numerically the second order corrections appear to be very moderate: 
Tt{nn,/mt = 0) = To (l - 0.866 a,(mt) - 1.69(5) af) . (18) 



For nib/nit = 4.8/175 we find 

a2{nib/mt = 4.8/175) = -12.4(3), a^^^{mb/nit = 4.8/175) = -16.3(2); 
Ttimb/mt = A.S/ 175) = 0.998 x Fq (l - 0.857a,(mf) - 1.68(4)«^) . (19) 

Using as{nit) = 0.11, we find that the 0{al) corrections to the top decay width are 
very close to 2%. 

Finally, we list here also the values obtained for 6 = 0.7, or nib/ nit = 0.3, which 
after replacement t,b ^ b,c are relevant for the differential width of 6 — > c + leptons, 
when the leptons are emitted parallel to each other, i.e. with a zero invariant mass: 
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For rrib/mt = 0.3 : 

ap = 1.37261(2), ua = -3.09498(4), 

ttL = 1.31091(5), an = -0.05061395(1). (20) 

At rrih/mt = 0.3 one has to account for the finite b mass in the virtual loop. This 
amounts to replacing an by an, which corresponds to S^Afj in At 5 = 

1 — rrib/mt = 0.7 we find 

aH = aH + (Woop) ~ -0.0506 + 0.6333 = 0.5827. (21) 



4 Conclusion 

We have presented a calculation of the 0{al) corrections to the decay width of the 
top quark T(t ^ bW) in the limit rrit ^ mw For asirrit) = 0.11, we found (see 
eqs. (0, ^)) that the second order QCD corrections decrease the value of the top 
width by 2%. 



To perform this calculation, we used the methods described in detail in [^. Also, 
since the process t — > bbbW contributes to the top decay width at order 0{al), we 
had to develop a new technique to deal with the "non-planar" interference diagrams. 

The decay rate of the top quark is proportional to the third power of its mass. In our 



calculations we used the pole mass of the top quark. It has been demonstrated |18] 
that the convergence of the perturbation series is improved if one parametrizes the 
width formula in terms of the MS mass. We have not performed this reparametriza- 
tion here, since the second order corrections are small even if the pole mass is used. 
We note, however, that employing the MS mass is likely to decrease them even 
further. 

With the electroweak corrections and effects of the W width on F^ known, our 
result provides the last missing ingredient in predicting the top decay width in the 
SM with accuracy ~ 1%. The remaining uncertainty is presently dominated by the 
top quark mass. It is difficult to say at the moment whether or not the experimental 
determination of the top quark width can be performed with a comparable precision. 
A future muon collider might be the best place for such measurements. 
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